TECHNICAL 

LIBRARY 


137536 

_  — .  .  O  l .  T-r-i  -i 


AD 


TECHNICAL  REPORT  ARBRL-TR-02544 


DERIVATION  OF  THE  P2T  CRITERION  (IN  CLOSED 
FORM)  BASED  ON  THE  FREY  SHEAR  BAND 

MECHANISM 


Evan  Harris  Walker 


January  1984 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER 

BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Approved  for  public  release;  distribution  unlimited. 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers  '  names  in  this  report 
does  not  constitute  indorsement  of  any  convnerdat  product. 


imilassjleIed 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

TECHNICAL  REPORT  ARBRL-TR-02544 

3.  RECIPIENT’S  CATALOG  NUMBER 

A.  TITLE  (and  Subtitle) 

DERIVATION  OF  THE  P2T  CRITERION  (IN  CLOSED  FORM) 
BASED  ON  THE  FREY  SHEAR  BAND  MECHANISM 

s.  type  of  report  a  period  covered 

Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR^ 

Evan  Harris  Walker 

8.  CONTRACT  OR  GRANT  NUMBER^ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory,  ARDC 

ATTN:  DRSMC-BLT (A) 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

AREA  &  WORK  UNIT  NUMBERS 

1L161102M43 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  AMCCOM,  ARDC 

Ballistic  Research  Laboratory,  ATTN:  DRSMC-BLA-S(A) 
Aberdeen  Proving  Ground,  MD  21005 

12.  REPORT  oate 

January  1984 

I  13.  HUMBER  OF  PAGES 

23 

1*.  MONITORING  AGENCY  NAME  &  ADDRESSfif  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (o t  thle  report) 

UNCLASSIFIED 

IS©.  DECLASSI  FI  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  rover#©  aide  If  neceeeary  and  Identify  by  block  number) 

Explosive  Sensitivity,  Critical  Energy,  Shear  bands  in  Explosives,  Detonation 
Criterion,  Composition-B,  Di-constituent  explosives 

20.  ABSTRACT  fCbntfiiu.  an  reveree  afdt  tt  neceeeery  amt  Identify  by  block  number ■)  (d  lC  )  2 

We  obtain  an  approximate  closed  form  equation  for  the  critical  P  T  criterion 
for  explosive  initiation  based  on  the  shear  band  theory  of  explosive  initiation. 
The  theory  assumes  that  shock  loading  of  explosives  produces  a  calculable  shear 
velocity  resulting  in  shear  band  formation.  Heating  within  the  shear  bands 
leads  to  hot  spot  formation  which  in  turn  leads  to  explosive  initiation.  An 
example  of  the  use  of  the  derived  equation  for  Composition-B  is  given. 

DO  »JAM*73  M73  EDITIOM  OF  »  MOV  65  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Deta  Entered) 


TABLE  OF  CONTENTS 


PACE 


I .  INTRODUCTION .  5 

II.  WORK  DONE  BY  SHEAR  FORCES .  7 

III.  THE  SHEAR  VELOCITY  IN  DI-CONSTITUENT  EXPLOSIVES  .  8 

IV.  THE  P2T  CONDITION . 11 

V.  CONCLUSIONS . 14 

REFERENCES . 15 

APPENDIX  A . 17 

DISTRIBUTION  LIST . 21 


3 


I .  INTRODUCTION 


In  a  previous  report/  an  expression  for  the  shear  band  shear  velocity  in 
a  shock  loaded  homogeneous  di-constituent  explosive  was  derived  and  used  with 
the  Frey  shear  band  theory  of  explosive  initiation  to  calculate  values  of 
2 

P  T  for  explosives,  where  P  is  the  pressure  to  which  the  explosive  is  shock 
loaded  and  T  is  the  temporal  duration  of  the  applied  pressure.  In  that 

2 

treatment,  the  origin  of  the  P  T  condition  as  such  was  not  derived.  The 
approach  used  a  computer  code  developed  by  R.  Frey2  (Explosive  Effects  Branch, 
Terminal  Ballistics  Division,  Ballistic  Research  Laboratory)  to  solve  the 
problem  of  heat  evolution  and  exchange  during  the  growth  of  shear  bands.  In 
the  present  paper,  approximate  expressions  are  obtained  for  the  development 
of  the  shear  bands  which  are  used  to  obtain  a  closed  form  expression  for  the 
2 

P  T  criterion.  Such  an  expression,  though  approximate,  has  great  utility  in 
exhibiting  the  important  parameters  that  determine  the  sensitivity  of  an  ex¬ 
plosive  in  terms  of  the  physical  and  thermodynamic  properties  of  that  explosive. 


2 

The  P  T  criterion  for  the  onset  of  detonation  was  first  proposed  by  Walker 
and  Wasley.3  The  criterion  has  been  extensively  studied  and  supported  by  ex¬ 
perimental  data.1*”9  The  criterion  is  essentially  equivalent  to  the  critical 
total  energy  fluence  £  (since  £  =  put;  u  being  the  particle  velocity,  and 

2 

p  =  P0uv,  pq  being  the  ambient  explosive  density,  we  obtain  E  =  p  t/pQv; 

p  v  is  nearly  constant  for  pressures  of  interest) .  An  understanding  of  the 

processes  by  which  initiation  of  detonation  occurs  requires  an  understanding 
2 

of  the  P  T  mechanism. 


-  .  2 

iE.  H.  Walker,  "Use  of  the  Shear  Band  Initiation  Model  to  Compute  Critical  P  T 
Values  for  Explosives BRL  Technical  Report  (in  preparation). 

2R.  B.  Frey,  "The  Initiation  of  Explosive  Charges  by  Rapid  Shear, "  Seventh 
Symposium  (International)  on  detonation,  pp  53-59,  16-19  Jun  81,  Annapolis,  MB. 

ZF.  E.  Walker  and  R.  J.  Wasley,  "Critical  Energy  for  the  Shock  Initiation  of 
Heterogeneous  Explosives,"  Explosives toffe  17,  9  (1969) . 
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7R.  C.  Weingart,  R.  K.  Jackson,  C.  A.  Honodel,  and  R.  S.  Lee,  "Shock  Initiation 
of  PBX9404  by  Electrically  Driven  Flyer  Plates,  "  Tenth  Symposium  on  Explosives 
and  Pyrotechnics,  preprint  UCRL-82162,  Jan  79. 

8y.  de  Longueville,  C.  Fauquignon,  and  H.  Moulard,  "Initiation  of  Several  Con¬ 
densed  Explosives  by  a  Given  Duration  Shock  Ware, "  Sixth  Symposium  (Inter¬ 
national)  on  Detonation,  pp  105-114,  ACR-221,  US  Government  Printing  Office, 


Washington,  D.C.,  20402.  ,  .  •  i 

9E.  H.  Walker,  "Computer  Model  of  Shield  Design  for  Artillery  Shell  Fratricide 

Suppression,"  BRL  Technical  Report  (in  preparation). 


A  number  of  studies*®  have  sought  to  derive  or  calculate  values  for  the 

2 

P  T  criterion,  but  none  has  provided  a  strictly  physical  theory  free  of  ad 
hoc  assumptions  as  to  the  origin  and  distribution  of  initiation  sites,  "hot 
spots,"  or  the  assumed  initiation  temperature  distribution  behind  the  shock 
front.  The  necessity  of  the  "hot  spot"  concept  of  Bowden1 5  and  Eyring1 6 
comes  from  the  fact  that, under  shock  initiation  conditions,  the  bulk  tempera¬ 
ture  of  the  explosive  produced  by  adiabatic  compression  is  generally  much 
too  low  to  ignite  the  explosive.  Bowden  and  Yoffee15  developed  the  theory 
of  adiabatic  compression  of  air  cavities  in  the  explosive  as  the  source  of 
the  hot  spots.  (See  also  the  work  of  Starkenberg.1 7  )  Mader1  8  has  suggested 
void  closure  in  the  explosive  as  the  mechanism  for  hot  spot  formation  and 
Seeley19  has  suggested  jetting  during  void  closure.  Delpuech,  et  al,  have 
proposed  shock  induced  free  radical  formation  as  the  mechanism  underlying 
hot  spot  formation.  Nevertheless,  these  mechanisms  have  not  been  tied  to  a 

P  T  criterion.  The  shear  band  mechanism  proposed  by  Frey,  however,  has 

been  used  to  calculate* 1 * * * * *  the  critical  value  of  P  T  for  Composition-B.  It  is 
therefore  important  to  tie  the  shear  band  mechanism  for  hot  spot  formation 

2 

directly  to  the  condition  that  the  critical  value  of  P  T  for  an  explosive  is 
a  constant . 


1 6  See  References  3,  6,  8. 
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PETN,"  Sixth  Symposium  (International)  on  Detonation,  pp  20-28,  ACR-221,  US 
Government  Printing  Office,  Washington,  D.C.,  20402. 
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Printing  Office,  Washington,  D.C.,  20402. 
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natior^pp  219-225,  16-19  Jun  81,  Annapolis,  MD. 

1  5F.  P.  Bowden  and  A.  D.  Ioffe,  Initiation  and  Growth  of  Explosions  in  Liquids 
and  Solids,  Cambridge  University  Press  (1952). 

16H.  Eyring,  R.  E.  Powell,  G.  H.  Duffrey,  and  R.  B.  Darlin,  Chemical  Review  45, 


69  (1949)  . 

i 7 j'  starkenberg j  u Ignition  of  Solid  High  Explosive  by  the  Rapid  Compression  of 

an  Adjacent  Gas  Layer,"  Seventh  Symposium  (International)  on  Detonation,  pp  1-12, 
16-19  Jun  81,  Annapolis,  MD. 

18C.  L.  Mader,  Phys  Fluids  8_,  10  (1965).  _  7 

19£.  B.  Seeley,  "A  Proposed  Mechanism  for  Shock  Initation  of  Low  Density  Granutar 

Explosives,"  Proceedings  of  the  Fourth  Electric  Initiator  Symposium  (1963). 

20 A.  Delpuech,  J.  Cherville,  and  C.  Michaud,  "Molecular  Electronic  Structure  and 
Initiation  of  Secondary  Explosives , "  Seventh  Symposium  (International)  on 

Detonation ,  pp  36-44,  16-19  Jun  81,  Annapolis,  MD. 
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II.  WORK  DONE  BY  SHEAR  FORCES 


The  work  W  per  unit  volume  done  in  shearing  is  given  by  the  general 
expression 

dW  =  a-  •  d  £•  •  Cl) 

Y  j 

where  a.,  is  the  shear  stress  in  the  plane  normal  to  the  x.  -  axis  directed 
ij  1 

along  the  x.  -  axis  of  a  Cartesian  coordinate  system,  and  e..  is  the  corre- 
J  . 

sponding  strain.  If  we  have  shears  only  on  surfaces  normal  to  the  y  -  axis, 
directed  only  in  the  z  direction,  with  the  corresponding  strains,  we  can  write 
simply 


dW 


0  d  e 
yz  yz 


(2) 


The  shear  stress  for  a  viscoplastic  material  for  such  a  flow  is 


+  0 

S 


(3) 


where  vz  is  the  shear  velocity  in  the  direction  of  the  z  -  axis  varying  in  the 

y  direction.  The  strength  of  the  material  is  given  by  as,  while  y  is  the 

viscosity.  The  development  of  the  shear  bands  takes  place  in  essentially  two 
stages:  (l)initial  development  of  shear  failure  planes  and  (2)  viscous  heating  of 
the  shear  layers.  In  the  first  stage,  the  material  undergoes  elastic  deformation 
up  to  shear  failure,  so  that  the  viscosity  term  in  Equation  (3)  plays  no  role; 
and  in  the  second  stage  the  shear  strength  drops  to  zero: 


(3) 


a=o  (4) 

s 

During  stage  two,  the  principal  heating  phase,  we  can  re-express  Equation 
as 


d^ 


,  ,,  6 (dZ/dt)  Sy 

CX„„  =  ydv^/dy  =  y  --  Y,  ■  •  ■  =  y 


yz 


Sy 


or  writing  e  for  Sz/Sy  and  rearranging 
yz 


de  =  —  a  dt 
yz  y  yz 


(5) 


(6) 


7 


Therefore,  we  can  substitute  for  de 

’  yz 


in  Equation  (2)  to  obtain,  on  integrating 


(7) 


III.  THE  SHEAR  VELOCITY  IN  DI -CONSTITUENT  EXPLOSIVES 


Walker1  has  given  an  expression  for  the  shear  velocity  in  an  explosive 
consisting  of  two  materials.  For  a  matrix  material  having  density  and  sound 

velocity  u^  containing  grains  with  density  p ^  and  sound  velocity  a  shock 

pressure  p  gives  rise  to  particle  velocities  v^  and  v^  satisfying 


P1U1 


V1  = 


p2  U2  V2 


(8) 


Therefore,  the  shear  velocity  v  ,  the  difference  between  v^  and  v^  is 


As  shown  in  Figure  1,  the  shear  velocity  initially  causes  a  strain  to 
develop  in  the  matrix  material.  (A  similar  process  occurs  in  the  grains.)  Where 
2n  is  the  number  of  shear  bands  in  the  spacing  of  2D  between  grains  in  the 
explosive  mixture,  at  a  time  t  the  strain  n6d  (as  shown  in  Figure  1)  will  be 


n6d  =  v  t 
s 


(10) 


and  also  the  number  of  bands  between  the  wall  and  center  between  grains  is 
given  in  terms  of  the  average  band  spacing  d  by 


n  =  D/d 


(ID 


When  the 


strain  6d/d  reaches  a  critical  distortion  angle  0 


max’ 


6d/d  =  0 


max 


(12) 


8 
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Figure  1  -  Schematic  diagram  showing  (a)  the  differential  velocity  field  produced  behind  a  shock  passing 
through  a  di-component  explosive,  (b)  a  detail  of  the  initial  strain  due  to  the  differential 
velocity  flow  (exaggerated) ,  and  (c)  a  schematic  of  the  shear  bands  of  thickness  x^  separated 
by  unsheared  material  of  thickness  d. 


shears  develop  in  the  material.  As  shears  begin  to  develop,  however,  stress 
relief  waves  propagate  from  these  incipient  shears  into  the  surrounding  material. 
This  means  that  the  next  nearest  shear  will  be  located  on  the  average  a  distance 
Ujt  away.  Since  this  distance  is  also  d  we  have 


u,  t  =  d 
1  » 


so  that  solving  we  have 


d  = 


0 

max 


DUj/v 


s 


(13) 


(14) 


0  can  be  re-expressed  in  terms  of  Poisson's  ratio  v,  the  ultimate  shear 
max  r  9 

strength  and  the  Young 's  elastic  modulus  E  by 


0 

max 


2(1  +  v) 


a  /E 
max 


(15) 


The  width  of  the  band  when  a  —>  0  can  be  calculated  in  terms  of  the  heat 

s 

made  available  by  the  formation  of  these  shears  during  the  time  t  in  which  the 

shear  banding  structure  forms.  Thus,' using  energy  conservation  to  determine 
the  thickness  of  the  layer  heated,  we  write. 


xbpiQf 


Cf  v,  t 
max  b 


(16) 


where  is  the  latent  heat  of  fusion  and  V^,  the  shear  velocity  for  a  single 
band,  is  given  by 


Thus,  we  obtain  for  x^ 


vb  =  vs/n 


a 

max 


n  p 


0 

max 


1  ^f 


D 


(17) 


(18) 


or,  substituting  from  Equation  (15)^ 

nxb  =  2  (l+v)Damax2/PlQfE 


(19) 


10 


From  Equations  (9)  and  (19) 

the  shear:  dy 

z  z 

3y"  ~  nx^  =  P 

This  expression  can  then  be 
tion  (4)  holds. 


we  can  obtain  the  average  velocity  gradient  across 

i  32U2  -  P1U1 

(20) 


2 (1+v) Da 


max 


P1U1P2U2 


used  to  obtain  O ^  from  Equation  (3)  where  Equa- 


IV.  THE  P  T  CONDITION 


Substituting  from  Equation  (20)  into  Equation  (3)  with  0$  set  to  zero, 
we  have 


o  =  p 
yz  r 


yp^Q^E  /  p2u2  “  Piui 


2 (1+v) Da 


max 


P1U1P2U2 


(21) 


Substituting  this  expression  for  a into  Equation  (7)  gives 


W  =  y 


plQfE 


2 (1+v) Da 


max 


P2U2  "  P1U1 

P1U1P2U2 


/ 


p2dt 


(22) 


where  the  constants  (or  slowly  varying  parameters)  have  been  factored  from 
under  the  integral . 

The  energy  W  (per  unit  volume)  goes  into  heating  the  material  in  the 
shear  band.  Where  At.  is  the  temperature  rise  above  the  melt  temperature 

required  to  bring  the  matrix  to  the  "ignition  temperature,"  we  can  express  the 
critical  value  for  W  by 


W 


crit 


p,C  AT. 

1  v  lg 


(23) 


where  C  is  the  specific  heat  of  the  matrix  material.  Strictly  speaking,  there 

is  no  specific  ignition  temperature.  Instead,  there  exists  a  relationship, 
the  Arrhenius  equation,  relating  rate  of  heat  evolution  to  the  temperature  o 
the  explosive.  However,  for  the  shock  loading  times  of  major  interest,  i.e., 
in  the  range  of  1-100  ys,  the  dependence  of  T.  on  the  shock  loading  time 

varies  slowly.  Heat  evolution  rises  exponentially  with  temperature,  which  in 
turn  further  elevates  the  explosive  temperature;  as  a  result,  the  initial  _ 
temperature  necessary  to  produce  complete  combustion  within  the  shear  band  in 
a  time  of  lys  will  differ  very  little  from  that  required  for  complete  combus¬ 
tion  in  100  ys.For  this  reason  we  can  ignore  the  dependence  of  Tig  on  shock 

loading  time  in  order  to  retain  a  closed  form  expression  for  Defining 


Qh 


as  the  P 


2 


time  integral  for  the  shear  band  heating  phase,  we  have  from 


Equations  (22)  and  (23) : 


11 


2 


2 


C  AT. 
v  ig 

PXV 


2(1+v)Dct 


max 


QfE 


P1U1P2U2 

P2U2"P1U1 


(24) 


Equation  (24)  provides  an  approximate  expression  for  the  P  T  criterion  in  terms 
of  the  physical  and  thermodynamic  properties  of  an  explosive  under  the  conditions 
of  matrix  initiation.  A  corresponding  expression  for  grain  initiation  may  be 

2 

given  in  which  D  and  amax  are  calculated  for  grain  shearing.  The  P  T  criterion 

is  given  by  the  expression  in  Equation  (24)  if  the  matrix  goes  into  shear  before 
the  grains,  as  in  Compos it ion- B.  If  the  grains  shear  with  the  matrix  deforming, 
the  quantities  p2,u2  will  be  for  the  matrix  and  all  other  quantities  for  the 
grains . 


Matrix  (TNT)  C 


U] 

y 

V 


1.29  x  105  J/Kg°K 


AT.  =  T. 
ig  ig 


T  .  »720°C 
melt 


1.63  x  103  Kg/m3 
2.08  x  103  m/s 
1.39  x  10  ^  Kg/ms 


=  0.3 


For  Composition-B  we  have  the  values: 

D  =  5  x  lO'V*15 


max 


7.0  x  10  Pa 


=  8.96  x  10  JKg 


E  =  5.4  x  10  Pa 


Grain  (RDX) 


p2  -  1.8  x  103  Kg/m3 

u  =  2.65  x  103m/s 


These  values  yield 
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9.8  x  10  Pas 


(25) 


In  addition  to  the  frictional  heating  in  the  shear  bands,  once  ignition 
of  the  explosive  at  these  hot  spots  has  occurred,  combustion  must  continue  if 
a  detonative  shock  pressure  is  to  form.  If  a  rarefaction  wave  arrives  too 
early,  in  the  process,  the  rate  of  heat  evolution  will  drop  and  cooling  or 
mechanical  dispersal  of  the  explosive  can  extinguish  the  ignition  process. 
Therefore,  if  the  burn  is  to  be  maintained,  the  externally  applied  shock 
pressure  must  be  maintained  until  this  external  source  of  pressure  is  re¬ 
placed  by  internal  sources.  The  internal  sources,  of  course,  must  be  the 
hot  spots  which  generate  pressure  by  combustion.  We  require  that  the  ex- 
%  email  y  applied  pressure  p  be  maintained  for  a  time  t^  sufficiently  long  that, 
^  See  Appendix  A. 
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if  at  the  end  of  this  time  the  external  pressure  is  relieved,  the  internal 
pressure  will  just  match  that  loss  of  pressure. (Note  that  this  internal  pres¬ 
sure  will  not  lead  to  mechanical  disruption  of  the  explosive  until  a  relief 
wave  can  arrive  from  the  explosive  surface.)  Since  we  are  concerned  with 
compact  explosives  here,  the  energy  density  is  equal  to  the  pressure  generated. 
Ihus,  to  produce  a  pressure  p  from  the  combustion  of  an  explosive  having  a 
"bomb"  pressure  (bomb  at  full  density)  p*,  a  fraction  p/p*  must  be  burned. 

Where  is  the  burn  velocity  and  d  the  distance  between  shear  bands,  we 
require 


p/p*  =  Ub  tb/(l/2  d)  (26) 

2 

Thus,  the  total  P  T  value  of  the  critical  condition  for  explosive  initiation 
will  be 


Qtot  =  %  +  %  (27) 

where  Qb  is  given  by 

Qb  "  p2  tb  =  P3(i/2Ub  P*  (28) 

Of  course,  the  burn  velocity  Ub  is  itself  a  function  of  the  pressure. 

Kondrikov ,  Raikova,  and  Samsonof22  give 

Ub  =  PBPV  (29) 

where  p  is  the  explosive  density,  p  the  pressure,  B  and  v  are  constants  with 

values  B=0.0219g/cm2s  atmV(20-100  Atm),  =  0.0094g/cm2s  atmV(100-450  Atm); 

V  =  0.77  (20-100  Atm) ,  -  0.95  (100-450  Atm).  We  obtain  p*  for  TNT  using 

the  TIGER  computer  code,23  giving  p*  =  9.67  x  109  Pa.  At  10  kbar, 

(2.03  x  108Pa),  Ub  =  0.21  m/s  so  that  Qb  =2.79  x  10U  Pa2s.  Since  Ub  is 

approximately  proportional  to  pressure,  Qb  varies  quadratically  with  the 

pressure.  Thus,  Q,  will  not  be  as  important  as  Q,  until  we  reach  shock 
d  h 

pressures  of  more  than  6.75  x  10  Pa  (67  kbar).  This  latter  value  is  only 

approximate  since  this  pressure  lies  outside  the  range  for  which  values  of 

2  1° 
are  accurate.  The  critical  P  T  value  obtained  experimentally  is  4.41  x  10 

2  2 
Pa  s.  Thus,  the  theoretical  value,  Q  ,  for  the  P  T  criterion  as  given  in 

Equation  (25)  is  reasonably  good. 

22 

B.  N.  Kondrikov,  V.  M.  Raikova,  and  B.  S.  Samsonof,  Fizika  Goreniya  i 
Vzryva,  9_  #1,  84-90  (1973). 

23 

M.  Cowperthuaite  and  W.  H.  Zwisler,  TIGER  Computer  Program,  SRI  Inter¬ 
national  Publication  No.  Z106,  SRI  International ,  333  Ravenswood  Ave, 

Menlo  Pk,  CA  94025. 
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V.  CONCLUSIONS 


We  have  obtained  an  approximate  closed  form  expression  for  the  P2T  cri¬ 
terion  for  explosive  initiation  based  on  the  shear  band  theory  of  explosive 
ignition.  While  the  resulting  expression  depends  on  quantities  that  are  only 
poorly  known  in  many  cases,  the  equation  exhibits  explicitly  the  factors  con¬ 
trolling  explosive  sensitivity.  The  resulting  equation  shows  a  dependence  of 
2 

the  P  T  criterion  on  particle  size  for  fixed  mass  fractions  of  the  constituents, 

2 

and  for  fixed  particle  size,  shows  P  T  increases  as  the  portion  of  matrix 

2 

(TNT  in  Comp  B)  in  increased .  The  P  T  criterion  equation  obtained  here 
exhibits  explicitly  the  important  parameters  determining  the  sensitivity  of 
an  explosive  in  terms  of  its  physical  and  thermodynamic  properties. 
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APPENDIX  A 


CALCULATION  OF  AVERAGE  GRAIN  SIZE  OF  RDX  IN  COMPOS ITION-B 


The  value  of  D,  the  average  grain  diameter  used  in  Equation  (25)  was  obtained 
from  data  for  the  composition-B  shown  in  Figure  2  of  Reference  1.  The  data  are 
given  in  Table  A1 . 


TABLE  A1 

Percent  of  RDX  having  Grain  Sizes  in  the  Indicated  Size  Ranges 
Percent  RDX  (%)  Size  Range  (ym) 


2.0 

> 

900 

9.8 

> 

350 

35.3 

> 

200 

39.7 

> 

100 

13.2 

< 

100 

To  use  these  data  to  obtain  D,  we  must  make  an  assumption  as  to  the  packing  of 
the  grains  and  their  shape.  We  assume  the  grains  to  be  spherical  and  equally 
spaced  in  a  close  packed  hexagonal  arrangement.  This  does  not  mean  that  the 
grains  are  in  contact  but  that  each  grain  of  radius  rg  at  the  center  of  a 
sphere  of  radius  R,  which  we  can  treat  as  pseudograins  in  contact  with  each 
other  in  the  close  packed  hexagonal  arrangement.  The  volume  of  a  grain  is 


V 

gram 


The  volume  of  a  pseudograin  is 


pseudo 


and  the  quantities  R,  r 


g’ 


and  D  are  related  by 


(Al) 


(A2) 


D  =  R  -  r  (A3) 

g 

The  volume  of  a  cell  containing  one  pseudograin  for  the  close  packed  hexagonal 
arrangement  (boundaries  of  the  cell  are  defined  by  the  tangent  planes  at  the 
points  of  contact  between  spheres)  is 


cell 


=  yr-r 


(A4) 
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Therefore,  g  is  given  by 


g  =  Vain/Vcell  ’  '  V/3  ^ 


(A5) 


For  the  mass  fraction,  we  note  that  the  masses  of  the  grain,  Mgr&in,  and  matrix, 

M  .  ,  are  given  in  terms  of  the  respective  densities,  Pgra^n  an<^  pmatrix’  ^ 
matrix  ® 


w  4  3 

M  =  —  -nrp  r 

grain  3  g  g 


(A6) 


and 


^matrix  ^matrix  ^  R  -  y  *  rg  J 


(A7) 


With  the  mass  fraction  of  the  grain  material  being  gm,  we  have 

3 


Sm 


w  P  v 
g  g 


vr 


matrix 


.  R  +  t  /( 
IX  I 


pgrain  p matrix)  Vg 


ix^ 


(A8) 


From  the  data  in  Table  A1  we  can  obtain  the  weighted  value  of  rgr&in  : 


=2>i 


r  r 
g  g 


(5.65  x  10"5m)  2 


X^w.  r 

i  g 


(A9) 


where  the  term  1/r  represents  the  grain  weights,  w.  their  percentages  and 

j.  ^  used  since  we  seek  the  value  of  D  for  use  in  Equation  (24)  .  Where 

gg  =  0.6,  p  .  =1.80  g/cm3,  and  p  .  =1.63  g/cm3.  Equations  (A8)  and 

&m  grain  matrix 

(A3)  can  be  solved  to  obtain 


D  =  4.94  x  10  ^m 


(A10) 


which  was  rounded  off  to  5  x  10  m  for  use  in  Equation  (24) 
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